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ABSTRACT 

If both the day-side and night-side effective temperatures of a planet can be measured, it is possible 
to estimate its Bond albedo, < Ab < 1, as well as its day-night heat redistribution efficiency, 
< e < 1. We attempt a statistical analysis of the albedo and redistribution efficiency for 24 
transiting exoplanets that have at least one published secondary eclipse. For each planet, we show 
how to calculate a sub-stellar equilibrium temperature, To, and associated uncertainty. We then use 
a simple model- independent technique to estimate a planet's effective temperature from planet/star 
flux ratios. We use thermal secondary eclipse measurements — those obtained at A > 0.8 micron — 
to estimate day-side effective temperatures, Td, and thermal phase variations — when available — to 
estimate night-side effective temperature. We strongly rule out the "null hypothesis" of a single Ab and 
e for all 24 planets. If we allow each planet to have different parameters, we find that low Bond albedos 
are favored (Ab < 0.35 at la confidence), which is an independent confirmation of the low albedos 
inferred from non-detection of reflected light. Our sample exhibits a wide variety of redistribution 
efficiencies. When normalized by To, the day-side effective temperatures of the 24 planets describe 
a uni-modal distribution. The two biggest outliers are GJ 436b (abnormally hot) and HD 80606b 
(abnormally cool), and these are the only eccentric planets in our sample. The dimensionless quantity 
Td/To exhibits no trend with the presence or absence of stratospheric inversions. There is also no 
clear trend between T^/Tq and To- That said, the 6 planets with the greatest sub-stellar equilibrium 
temperatures (T > 2400 K) have low e, as opposed to the 18 cooler planets, which show a variety 
of recirculation efficiencies. This hints that the very hottest transiting giant planets are qualitatively 
different from the merely hot Jupiters. We propose an explanation of this trend based on how a 
planet's radiative and advective times scale with temperature: both timescales are expected to be 
shorter for hotter planets, but the temperature-dependance of the radiative timescale is stronger, 
leading to decreased heat recirculation efficiency 

Subject headings: methods: data analysis — (stars:) planetary systems — 



1. INTRODUCTION 

Short-period exoplanets are expected to have atmo- 
spheric compositions and dynamics that differ signifi- 
cantly from Solar System giant planetfl These planets 
orbit ~ 100 x closer to their host stars than Jupiter does 
from the Sun. As a result, they receive ~ 10 4 x more flux 
and experience tidal forces ~ 10 6 x stronger than Jupiter. 
In contrast to Jupiter, which releases roughly as much 
power in its interior as it receives from the Sun, short- 
period exoplanets have power budgets dictated by the 
flux they receive from their host stars. Roughly speak- 
ing, the stellar flux incident on a planet does one of two 
things: it is reflected back into space, or advected else- 
where on the planet and re-radiated at different wave- 
lengths. The physical parameters that describe these 
processes are the planet's Bond albedo and redistribu- 
tion efficiency. 

1.1. Albedo 
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Giant planets in the Solar System have albedos greater 
than 50% because of the presence of condensed molecules 
(H2O, CH4, NH3, etc.) in their atmospheres. Planets 
with effective temperatures exceeding ~40 K should be 
cloud free, leading to albedos of 0.05-0.4 (jMarlev et al.1 
[1991 . If pressure-broadened Na and K opacity is impor- 
tant at optica l wave lengths (as it is for brown dwarfs, 
iBurrows et al.l 120001) . then the Bond albedos of ho t 
Jupiters may be less than 10% (Sud arsky et al.l [20001) . 
But the very hottest planets, the so-called class V extra- 
solar giant planets (T e g > 1500 K), might have very high 
albed os due to a high silicate cloud layer (Sudarsk y et al.l 
2000). For a planet whose albedo is dominated by 
clouds (as opposed to Rayleigh scattering) the albedo 
depends on the co mposition and size of cloud particles 
(|Seager et al.ll2000l ). 

Ear ly attempts to observe refl ec ted light from exoplan- 
ets 
[2002 



(ICharbonneau et al.| Il999t ICollier Cameron et all 
"at ILeigh et al.ll2003al lrl TRodrer et al.H2008l 12010ft b"T 



dicated that they might not be as reflective as Solar Sys - 
tem gas giants (for a review, see lLangford et all l2010h . 
Measurements of HD 209458b taken with the Cana- 
dian MOST sate llite revealed a very low albedo (< 8%, 
IRowe et al.ll2008ft , and it has since been taken for granted 
that all short-period planets have albedos on par with 
that of charcoal. 

From the standpoint of the planet's climate, the im- 
portant factor is not the albedo at any one wavelength, 
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A\, but rather the integrated albedo, weighted by the in- 
cident stellar spectrum, known as the Bond albedo and 
denoted in this paper as Ab- The relation between A\ 
and the planet's Bond albedo is not trivial. If the albedo 
is dominated by gray clouds, then the albedo at a sin- 
gle wavelength can indeed be extrapolated to obtain Ab- 
For non-gray reflectance spectra, however, it is critical to 
measure A\ at the peak emitting wavelength of the host 
star to obtain a good estimate of th e planet's energy bud - 
get. For example, as pointed out in lMarlev et al.l (1999). 
planets with identical albedo spectra, A\, may have radi- 
cally different Ab depending on the spectral type of their 
host stars. 

1.2. Redistribution Efficiency 

The first few measurements of hot Jupiter phase vari- 
ations showed signs that these planets are not a ll cut 
from the same cloth . lHarrington et al] (|2006l ) and 
iKnutson et al.l ()2007aD quoted very different phase func- 
tion amplitudes for the v Andromeda and HD 189733 
systems. It was not clear whether the differences were 
intrinsic to the planets, however, because the data 
were taken with different instruments, at different wave- 
lengths, and with very different observation schemes (in 
any case, subsequent re-analysis of the original data and 
newly aquired Spitzer observations of v Andromeda b 
paint a completely different picture of that system: 
ICrossfield et aIll201Qh . 

The uniform study presented in ICowan et al.l (|2007j ). 
on the other hand, showed that HD 179949b and 
HD 209458b exhibit significantly different degrees of heat 
recirculation, confirming suspicions. But it was not clear 
whether hot exoplanets were uni-modal or bi-modal in 
redistribution: are HD 179949b and HD 209458b end- 
members of a single distribution, or prototypes for two 
fundamentally different sorts of exoplanets? 

The presenc e or lack of a stratospheric tempera- 
ture inversion (iHubenv et al|[200l iFortnev et all 12001: 
iBurrows et ail 120071 l2008t iZahnle et al.l I2009D on the 
day-sides of exoplanets has been invoked to explain 
a purported bi-modality in recir culation efficiency on 
hot Jupiters ()Fortnev et al.ll2008D . The argument, sim- 
ply put, is that optical absorbers high in the atmo- 
sphere of extremely hot Jupiters (equilibrium temper- 
atures greater than ~ 1700 K) would absorb incident 
photons where the radiative timescales are short, mak- 
ing it difficult for these planets to recirculate energy. The 
most robust d etection of this temp erature inversion is for 
HD 209458b (|Knutson et al.l 120081 ). but this planet does 
not exhibit a large d ay-night brightness contrast at 8 /zm 
(|Cowan et al.l [20071. So while temperature inversions 
seem to exist in the majority of hot Jupiter atmospheres 
(jKnutson et al.|[20Tol) , their connection to circulation ef- 
ficiency — if any — is not clear. 

1.3. Outline of Paper 

It has been sugg ested (e.g., lHarrington et al.l 120061 : 
ICowan et al.l 120071) that observations of secondary 
eclipses and phase variations each constrain a combina- 
tion of a planet's Bond albedo and circulation efficiency. 
But observations — even phase variations — at a single 
waveband do little to constrain a planet's energy bud- 
get. In this work we show how observations in different 



wavebands and for different planets can be meaningfully 
combined to estimate these planetary parameters. 

In § 2 we introduce a simple model to quantify the 
day-side and night-side energy budget of a short-period 
planet, and show how a planet's Bond albedo, Ab, and 
redistribution efficiency, e, can be constrained by ob- 
servations. In § 3 we use published observations of 
24 transiting planets to estimate day-side and — where 
appropriate — night-side effective temperatures. We con- 
struct a two-dimensional distribution function in Ab and 
e in § 4. We state our conclusions in § 5. 

2. PARAMETERIZING THE ENERGY BUDGET 
2.1. Incident Flux 

Short-period planets have a power budget entirely dic- 
tated by the flux they receive from their host star, 
which dwarfs tid al heatin g or r emnant heat of forma- 
tion. Following Hansen (200H), we define the equi- 
librium temperature at the planet's sub-stellar point: 
T (t) = T cS (^/r(t)) 1 / 2 , where T cS and i?* are the star's 
effective temperature and radius, and r(t) is the planet- 
star distance (for a circular orbit r is simply equal to the 
semi-major axis, a). For shorthand, we define the geo- 
metrical factor a* = a /R*, which is directly constrain ed 
by transit lightcurves (jSeager fc Mallen-Orne las 2003) . 

The incident flux on the planet is given by Fi nc = 
JjCtbTq , and it is significant that this quantity has some 
associated uncertainty. For a planet on a circular orbit, 
the uncertainty in Tq — T s/^/a^ is related — to first 
order — to the uncertainties in the host star's effective 
temperature, and the geometrical factor: 
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(1) 



For a planet with non-zero eccentricity, To varies with 
time, but we are only concerned with its value at su- 
perior conjunction: secondary eclipse occurs at superior 
conjunction, when we are seeing the planet's day-side. 
At that point in the orbit, the planet-star distance is 
r sc = a(l — e 2 )/(l — esinw), where e and oj are the 
planet's orbital eccentricity and argument of periastron, 
respectively. 

For planets with non-zero eccentricity, the uncertainty 
in Tq is given by 



T 2 "^ST + 4a i \ !- 



( 



1 -p^ 



2(l-esinu) J esinui 



(2) 



where a e cos w and a e sin w are the observational uncertain- 
ties in the two components of the planet's eccentricitjQ. 

2.2. Emergent Flux 

At secondary eclipse, and in the absence of albedo or 
energy circulation, the equilibrium temperature of a re- 
gion on the planet depends on the normalized projected 

4 This formulation is preferable to an error estimate based on a e 
and (Tu, because the eccentricity and argument of periastron are 
highly correlated in orbital fits. That said, the uncertainties a e cos w 
and cr e s j n w are often not included in the literature, in which case 
we use a slightly different — and more conservative — formulation 
of the error budget using a e and <r^ . 
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distance, 7, from the center of the planetary disc as 
T(-/) = Xo(l — 7 2 ) 1 / 8 . The thermal secondary eclipse 
depth in this limit is given by: 

F^ = (Rp\ 2 ( Jw_\ 8 ( ho/xtai _ {] 
F, \Rj \XkT J \ J 

AXkTo/hcf dx 

X Jo exp^-VS)-!' (3) 

where T & * is the brightness temperature of the star at 
wavelength A. 

In the no-circulation limit, then, the day-side emer- 
gent spectrum is not exactly that of a blackbody, even 
if each annulus has a blackbody spectrum. But these 
differences are not important for the present study, since 
we are concerned with bolometric flux. By integrating 
Equation [3] over A, one obtains the effective tempera- 
ture of the day-side in the no-albedo, no-circulation limit: 
T E=n = (2/3) x / 4 T (see also lBurrows et al.ll200l iHansenl 
2008). Indeed, treating the planet's day-side as a uniform 
hemisphere emitting at this temperature gives nearly the 
same wavelength dependence as the more complex Equa- 
tion [3l The T e= o temperatures for our sample of 24 tran- 
siting planets are shown in Table 1. These set the max- 
imum possible day-side effective temperature we should 
expect to measure. 

The integrated day-side flux in the general — non-zero 
circulation — case is more subtle: heat may be trans- 
ported to the planet's night-side, and/or to its poles. In 
this paper we neglect the E-W asymetry in the planet's 
temperature map due to zonal flows and hence phase 
offsets in the thermal phase variations. Under this as- 
sumption, the day-night temperature contrast can more 
directly be extracted from the observed thermal phase 
variations. 

In practice, many studies have adopted a single param- 
eter to represent both zonal and meridional transport. It 
is instructive to consider the apparent day-side effective 
temperatures in various limits: uniform day-side temper- 
ature and T = on the night-side (this is often referred 
to as the planet's "equilibrium temperature"): T cqu = 
(1/2) 1 / 4 T ; in the case of perfect longitudinal transport 
but no latitudinal transport: Ti ong = (8/(37r 2 )) 1//4 T ; 
and in the limit of a uniform temperature everywhere 
on the planet: T uni = (l/£) l / A T . 

Comparing the apparent day-side temperatures in the 
three limits of circulation above leads to the following 
simple parametrization of the day-side effective temper- 
ature in terms of the planetary albedo, As, and circula- 
tion efficiency, e: 

(1 \ \ 1/4 
T d = T (l - A B )^ (± - A £ j , (4 ) 

where < e < 1. Note tha t e is related to — b ut dif- 
ferent from — the e used in (jCowan fc Agol|[2010l) . The 
former is merely a parametrization of the observed disk- 
integrated effective temperature, while the latter, which 
can take values from to 00, is a precisely defined ratio 
of radiative and advective timescales. The e = case is 
precisely equal to the e = case, while the e — > 00 limit 
is equivalent to e w 0.95. 
Our definition of e is similar to the iBurrows et al.l 



(2006) definition of P n and yields the same no-circulation 
limit. But our e = 1 limit produces a lower day-side 
brightness than the P n = 0.5 limit, because we as- 
sume that the planet's day-side has a uniform tempera- 
ture distribution in that limit (for a discussion of differ- 
ent redistribution param eterizations, see the appendix of 
ISpieeel fc BurrowsllMl . 

In reality, efficient longitudinal transport (read: fast 
zonal winds) may lead to more banding and therefore less 
efficient latitudinal transport. So one could argue that 
in the limit of perfect day-night temperature homoge- 
nization, both the day and night apparent temperatures 
should be T d = (8/(37r 2 )) 1 / 4 T , in between the Burrows 
et al. value of Ti = (l/3y/ 4 T and that suggested by 
our parameterization, T<j = (1/4) 1 / 4 T). At moderate 
day-night recirculation efficiencies, however, there is a 
good deal of latitudinal transport (I. Dobbs-Dixon, priv. 
comm.), so implicitly assuming a constant T oc cos 1 / 4 
latitudinal dependence — as done by Burrows et al. — is 
not founded, either. The bottom line is that any single- 
parameter implementation of advection is incapable of 
capturing the real complexities involved, but longitudi- 
nal transport is the dominant factor in determining day 
and night effective temperatures. 

Not withstanding the subtleties discussed above and 
noting that cooling ten ds to latitudinaly h omogenize 
night-side temperatures (| Cowan fc Ago1ll20lflh . we get a 
night-side temperature of: 

T n = T (l - A B ) 1/4 • (5) 

Note that T<j and T n are the equator- weighted tempera- 
tures of their respective hemispheres (ie, as seen by an 
edge-on viewer). As such, they will tend to be slightly 
higher than the hemisphere-averaged temperature, ex- 
cept in the e = 1 limit. This is also why the quantity 
Tj 4 + T 4 is still a weak function of e. 
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Fig. 1. — Different kinds of idealized observations constrain the 
Bond albedo, Ag and circulation efficiency, e, differently. A mea- 
surement of the secondary eclipse depth at optical wavelengths is 
a measure of albedo (solid line). A secondary eclipse depth at 
thermal wavelengths gives a joint constraint on albedo and recir- 
culation (dotted line). A measurement of the night-side effective 
temperature from thermal phase variations yields a constraint (the 
dashed line) nearly orthogonal to the day-side measurement. 

In Figure [1] we show how different kinds of observa- 
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tions constrain Ab and e. For this example, we chose 
constraints consistent with Ab — 0.2 and e — 0.3. The 
solid line is a locus of constant Ab] the dotted line is 
the locus of constant Td/T ; the dashed line is a lo- 
cus of constant T n /T . From this figure it is clear that 
the measurements complement each other: measuring 
two of the three quantities (Bond albedo, effective day- 
side or night-side temperatures) uniquely determines the 
planet's albedo and circulation efficiency. When obser- 
vations have some associated uncertainty, they define a 
swath through the Ab~s plane. 

3. ANALYSIS 
3.1. Planetary & Stellar Data 

We begin by considering all the photometric obser- 
vations of short-period exoplanets published through 
November 2010, summarized in Table [TJ We have dis- 
carded photometric observations of non-transiting plan- 
ets because of their unknown radius and orbital inclina- 
tion^]. This leaves us with 24 transiting exoplanets for 
which there are observations in at least one waveband 
at superior conjunction, and in some cases in multiple 
wavebands and at multiple planetary phases. 

Stellar and planetary data are taken from the Ex- 
oplanet Encyclopedia (exoplanet.eu), and references 
therein. We repeated parts of the analysis with the 
Exoplanet Data Explorer database (exoplanets.org) and 
found identical results, within the uncertainties. When 
the stellar data are not available, we have assumed typi- 
cal parameters for the appropriate spectral class, and so- 
lar metallicity. Insofar as we are only concerned with the 
broadband brightnesses of the stars, our results should 
not depend sensitively on the input stellar parameters. 

Knowing the stars' T e g, logg and [Fe/H], we 
use the PHOENIX/N extGen stellar spectrum grids 
(|Hauschildt et all |l999) to determine their brightness 
temperatures at the observed bandpasses. At each wave- 
band for which eclipse or phase observations have been 
obtained, we determine the ratio of the stellar flux to the 
blackbody flux at that grid star's T e g. We then apply this 
factor to the T e g of the observed star. 

It is worth noting that the choice of stellar model leads 
to systematic uncertainties in the planetary brightness 
that are of order the photometric uncertainties. For 
example, iChristiansen et al.1 (120101 ) use stellar models 
fo r HAT-P-7 from iKurucd (|2005D while we use those 
of lHauschildt et al.l (|1999tT The resulting 8 bright- 
ness temperatures for HAT-P-7b differ by as much as 
60 K, or slightly more than la. Our uniform use 
of lHauschildt et all (|1999l ) models should alleviate this 
problem, however. 

3.2. From Flux Ratios to Effective Temperature 

The planet's albedo and recirculation efficiency gov- 
ern its effective day-side and night-side temperatures, T<j 
and T n , respectively. Observationally, we can only mea- 
sure the brightness temperature, ideally at a number of 
different wavelengths: 71, (A). If one knew, a priori, the 



5 For completeness, these are: r-Bootis b, v- Andromeda b, 
51 Peg b, Gl 876d H P 75 289 b HP 179949b and H P 46375b 
llCharbonneau et al.l [19991; ICollier Cameron et al.l I2002bj ; 
Leigh et al. 2003a b; Harrington et al. 2006; Cowan et al. 2007; 
Scagcr & Pcming 2009; Crossficld et al. 2010; Gaulmc et al. 2010) 



emergent spectrum of a planet, one could trivially con- 
vert a single brightness temperature to an effective tem- 
perature. Alternatively, if observations were obtained at 
a number of wavelengths bracketing the planet's black- 
body peak, it would be possible to estimate the planet's 
bolometric flux and hence its effective temp erature in a 
mo del- independent way fe.g.. iBarma r] [200cl . 

We adopt the latter empirical approach of converting 
observed flux ratios into brightness temperatures, then 
using these to estimate the planet's effective tempera- 
ture. The secondary eclipse depth in some waveband di- 
vided by the transit depth is a direct measure of the ratio 
of the planet's day-side intensity to the star's intensity 
at that wavelength, ip(\) . Knowing the star's brightness 
temperature at a given wavelength, it is possible to com- 
pute the apparent brightness temperature of the planet's 
day side: 



Tb(A) 



he 
Xk 



log 1 



Jic/\kT*(\) 



(6) 



On the Rayleigh- Jeans tail, the fractional uncertainty 
in the brightness temperature is roughly equal to the 
fractional uncertainty in the eclipse depth; on the Wien 
tail, the fractional error on brightness temperature can 
be smaller because the flux is very sensitive to tempera- 
ture. 

By the same token, a secondary eclipse depth and 
phase variation amplitude at a given wavelength can be 
combined to obtain a measure of the planet's night-side 
brightness temperature at that waveband. 

Since the albedo and recirculation efficiency of the 
planet are not known ahead of time, it is not immedi- 
ately obvious which wavelengths are sensitive to reflected 
light and which are dominated by thermal emission. For 
each planet, we compute the expected blackbody peak if 
the planet has no albedo and no recirculation of energy, 
A £= o = 2898/T e= o /im. Insofar as real planets will have 
non-zero albedo and non-zero recirculation, the day side 
should never reach T e= q, and the actual spectral energy 
distribution will peak at slightly longer wavelengths. The 
coolest planet in our sample, Gl 436b, would exhibit a 
blackbody peak at A e= o = 3.1 /im, while the hottest 
planet we consider, WASP-12b, has A e= o = 0.9 /xm. 
In practice this means that ground-based near-IR and 
space-based mid-IR (e.g., Spitzer) observations are as- 
sumed to measure thermal emission, while space-based 
optical observations (MOST, CoRoT, Kepler) may be 
contaminated by reflected starlight. 

In Figure [H we demonstrate two alternative techniques 
to convert an array of brightness temperatures, lb (A), 
into an estimate of a planet's effective temperature, T e g. 
The solid black lin e shows a mod e l spe ctrum of ther- 
mal emission from iFortnev et al.l (|2008f ). with an ef- 
fective temperature of T c g = 1941 K shown with the 
black dashed line. The expected blackbody peak of 
the planet is marked with a vertical dotted line. The 
red points are the expected brightness temperatures in 
the J, H, and K s bands (crosses), as well as the IRAC 
asterisks) and MIPS (diamond) instruments on Spitze r 
Fazio et al.|[200l iRieke et al.|[200l iWerner et afll200l . 
Since the majority of the observations of exoplanets have 
been obtained with Spitzer IRAC, we focus on estimat- 
ing T e g based only on brightness temperatures in those 
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Fig. 2. — The sol id black line shows a model spectrum from 
Fortncv ct al. (2008) including only thermal emission (ic: no re- 
flected light). The planet's effective temperature is shown with the 
black dashed line, while the expected wavelength of the blackbody 
peak of the planet is marked with a black dotted line. The red 
points show the expected brightness temperatures in the J, H, and 
K s bands (crosses), as well as the IRAC (asterisks) and MIPS (dia- 
mond) instruments on Spitzer. The linear interpolation technique 
described in the text is shown with the red line. 

four bandpasses. 

Wien Displacement: The first approach is to simply 
adopt the brightness temperature of the bandpass clos- 
est to the planet's blackbody peak (the black dotted 
line). If only the four IRAC channels are available, the 
best one can do is the 3.6 /xm measurement, yielding 
T e g — 1925 K. There is — however — some subtlety in 
estimating the peak wavelength, as this is dependent on 
knowing the planet's temperature (and hence As and e) 
a priori. 

Linear Interpolation: The linear interpolation tech- 
nique, shown with the red line in Figure [2j obviates the 
need for an estimate of the planet's temperature. The 
brightness temperature is assumed to be constant short- 
ward of the shortest-A observation, and longward of the 
longest-A observation. Between bandpasses, the bright- 
ness temperature changes linearly with A. As long as 
the various brightness temperatures do not differ grossly 
from one another, this technique implicitly gives more 
weight to observations near the hypothetical blackbody 
peak. The bolometric flux of this "model" spectrum is 
then computed, and admits a single effective tempera- 
ture, which is T e ff = 1927 K for the current example. 
Since we hope to apply our routine to planets with well 
sampled blackbody peaks, we adopt the linear interpola- 
tion technique, as it can make use of multiple brightness 
temperature estimates near the peak. 

The two techniques described above produce similar ef- 
fective temperatures, though — unsurprisingly — neither 
gives precisely the correct answer. But these system- 
atic errors are comparable or smaller than the photo- 
metric uncertainty in observations of individual bright- 
ness temperatures (see Table 1). The best IR observa- 
tions for the nearest, brightest planetary systems (e.g., 
HD 189733b and HD 209458b) lead to observational un- 
certainties of approximately 50 K in brightness temper- 
ature. For many planets, the uncertainty is 100-200 K. 
By that metric, either the Wien displacement or the lin- 
ear interpolation routines give adequate estimates of the 
effective temperature, with errors of 16 K and 14 K, re- 



spectively. 

We make a more quantitative analysis of the systematic 
uncertainties involved in the Linear Interpolation tem- 
perature estimates as follows. We produce 8800 mock 
data sets: 100 realizations for 11 models and data in 
up to 8 wavebands (J, H, K, IRAC, MIPS; Since this nu- 
merical experiment chooses random bands from the eight 
available, the results should not be very different if ad- 
ditional wavebands are considered). We run our Linear 
Interpolation technique on each of these and plot in Fig- 
ure [3] the estimated day-side temperature normalized by 
the actual model effective temperature versus the num- 
ber of wavebands used in the estimate. The temperature 
estimates cluster near T est /T e g = 1, indicating that the 
technique is not significantly biased. The scatter in es- 
timates decreases as more wavebands are used, from a 
standard deviation of 7.6% if only a single brightness 
temperature is used, down to 2.4% if photometry is ac- 
quired in eight bands. We incorporate this systematic 
error into our analysis by adding it in quadrature to 
the observational uncertainties described in the follow- 
ing paragraph. This has the desirable effect that planets 
with fewer observations have a larger systematic uncer- 
tainty on their effective temperature. 
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Fig. 3. — The Linear Interpolation technique for estimating day- 
side effective as tested on a suite of eleven hot Jupiter spectral 
models provided by J.J. Fortney. The y-axis shows the estimated 
day-side effective temperature normalized by the actual model ef- 
fective temperature. The x-axis represents the number of bright- 
ness temperatures used in the estimate. Each color corresponds to 
one of the eleven models used in the comparison. The black error 
bars represent the standard deviation in the normalized tempera- 
ture estimates. 

In practice, we would like to propagate the photometric 
uncertainties to the estimate of T e s. For the Wien Dis- 
placement technique, this uncertainty propagates triv- 
ially to the effective temperature. For the linear inter- 
polation technique, a Monte Carlo can be used to esti- 
mate the uncertainty in T e g: the input eclipse depths 
are randomly shifted 1000 times in a manner consistent 
with their photometric uncertainties — assuming Gaus- 
sian errors — and the effective temperature is recomputed 
repeatedly. The scatter in the resulting values of T e g pro- 
vides an estimate of the observational uncertainty in the 
parameter, to which we add in quadrature the estimate 
of systematic error described above. The resulting uncer- 
tainties are listed in Table 1. These uncertainties should 
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be compared to the uncertainties in T £= q (also listed in 
Table 1), which are computed using the uncertainty in 
the star's properties and the planet's orbit. 

There are two practical issues with the linear interpo- 
lation temperature estimation technique. In some cases, 
only upper limits have been obtained, therefore one could 
set -0 = 0, with the appropriate 1-sigma uncertainty. But 
this approach leads to huge uncertainties in T c ff for plan- 
ets with a secondary eclipse upper-limit near their black- 
body peak. Instead of "punishing" these planets, we opt 
to not use upper-limits (though for completeness we in- 
clude them in Table 1). Secondly, when multiple mea- 
surements of an eclipse depth have been published for 
a given waveband, we use the most recent observation, 
indicated with a superscript "e" in Table 1. In all cases 
these observations either explicitly agree with their older 
counterpart, or agree with the re-analyzed older data. 

4. RESULTS 
4.1. Looking for Reflected Light 

For each planet, we use thermal observations (essen- 
tially those in the J, H, K s , and Spitzer bands) to es- 
timate the planet's effective day-side temperature, Td, 
and — when phase variations are available — T n . These 
values are listed in Table 1. In five cases (CoRoT- 
lb, CoRoT-2b, HAT-P-7b, HD 209458b, TrES-2b), sec- 
ondary eclipses and/or phase variations have been ob- 
tained at optical wavelengths. Such observations have 
the potential to directly constrain the albedo of these 
planets. One approach is to adopt the Td from thermal 
observations and calculate the expected contrast ratio at 
optical wavelength s, under the assumption of blackbody 
emission (see also iKipping fc Bakosll20Tol ). Insofar as 
the observed eclipse depths are deeper than this calcu- 
lated depth, one can invoke the contribution of reflected 
light and compute a geometric albedo, A g . If one treats 
the planet as a uniform Lambert sphere, the geometric 
albedo is related to the spherical albedo at that wave- 
length by A\ = \A g . These values are listed in Table 1. 

But reflected light is not the only explanation for an 
unexpectedly deep optical eclipse. Alternatively, the 
emissivity of the planets may simply be greater at op- 
tical wavelengths than at mid-IR wavelengths, in agree- 
ment with realistic spectral models of hot Jupiters, which 
predict brightness temperatures greater than T c g on the 
Wien tail (see, for example, the Fortney et al. model 
shown in Figure[2j which does not include reflected light). 
Note that this increase in emissivity should occur regard- 
less of whether or not the planet has a stratosphere: by 
definition, the depth at which the optical thermal emis- 
sion is emitted is the depth at which incident starlight 
is absorbed, which will necessarily be a hot layer - 
assuming the incident stellar spectrum peaks in the op- 
tical. 

Determining the albedo directly (ie: by observing re- 
flected light) can be difficult for short period planets, 
because there is no way to distinguish between reflected 
and re-radiated photons. The blackbody peaks of the 
star and planet often differ by less than a micron. There- 
fore, unlike Solar System planets, these worlds do not 
exhibit a minimum in their spectral energy distribution 
between the reflected and thermal peaks. The hottest 
— and therefore most ambiguous case — of the five tran- 



siting planets with optical constraints is HAT-P-7b. If 
one takes the mid-IR eclipse depths at face value, the 
planet has a day-side effective temperature of ~ 2000 K. 
When combined with the Kepler observations, one com- 
putes an albedo of greater than 50%. The large day-night 
amplitude seen in the Kepler bandpass is then simply 
due to the fact that the planet's night-side reflects no 
starlight, and the cool day-side can be attributed to high 
Ab and/or e. If, on the other hand, one takes the op- 
tical flux to be entirely thermal in origin (A\ = 0), the 
day-side effective temperature is ~ 2800 K. This is very 
close to that planet's T e= o, leaving very little power left 
for the night-side, again explaining the large day-night 
contrast observed by Kepler. The truth probably lies 
somewhere between these two extremes, but in any case 
this degeneracy will be neatly broken with Warm Spitzer 
observations: the two scenarios outlined above will lead 
to small and large thermal phase variations, respectively. 
It is telling that the only optical measurement in Table 1 
that is unanimously considered to constrain albedo - 
and not therm al emission— is th e MOST observations 
of HD 209458b (|Rowe et al.ll2008l ). the coolest of the five 
transiting planets with optical photometric constraints. 

The bottom line is that extracting a constraint on re- 
flected light from optical measurements of hot Jupiters is 
best done with a detailed spectral model. But even when 
reflected light can be directly constrained, converting this 
constraint on A\ into a constraint on ylg also requires 
detailed knowledge of both the star and the planet's spec- 
tral energy distributions, making for a model-dependent 
exercise. 

4.2. Populating the Aq-e Plane 

Setting aside optical eclipses and direct measurements 
of albedo, we may use the rich near- and mid-IR data to 
constrain the Bond albedo and redistribution efficiency 
of short-period giant planets. We define a 20 x 20 grid in 
Ab and e and use Equations 4 & 5 to calculate the nor- 
malized day-side and night-side effective temperatures, 
Td/To and T n /To, at each grid point, For each 

planet, we have an observational estimate of the day-side 
effective temperature, and in three cases we also have an 
estimate of the night-side effective temperature (as well 
as associated uncertainties). 

We first verify whether or not the observations are con- 
sistent with a single Ab and e. To evaluate this "null 
hypothesis", we compute the usual \ 2 = Si=i( mocie l — 
data) 2 /error 2 at each grid point. We use only the esti- 
mates of day-side and (when available) night-side effec- 
tive temperatures to calculate the x 2 , giving us 27-2=25 
degrees of freedom. The "best-fit" has x 2 = 132 (reduced 
X 2 = 5.3), so the current observations strongly rule out 
a single Bond albedo and redistribution efficiency for all 
24 planets. 

For 21 of the 24 planets considered here, we construct 
a two-dimensional distribution function for each planet 
as follows: 

PDF(i,j) = _^ =e -(T d -T d (:U)) 2 /(2^)\ (7) 
V 27rfJ I 

This defines a swath through parameter space with the 
same shape as the dotted line in Figure 1. 
For the three remaining planets (HD 149026b, 
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HD 189733b, HD 209458b), phase variation measure- 
ments help break the degeneracy: 




Albedo 

Fig. 4. — The global distribution function for short-period exo- 
planets in the Ag-s plane. The gray-scale shows the sum of the 
normalized probability distribution function for the 24 planets in 
our sample. The data mostly consist of infrared day-side fluxes, 
leading to the dominant degeneracy (see first the dotted line in 
Figure 1). 

We create a two-dimensional normalized probability 
distribution function (PDF) for each planet, then add 
these together to create the global PDF shown in Fig- 
ure 2J This is a democratic way of representing the data, 
since each planet's distribution contributes equally. 

In Figures [5] and [6] we show the distribution functions 
for the albedo and circulation of the 24 planets in our 
sample, obtained by marginalizing the global PDF of Fig- 
ure @] over either Ab or e. 
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Circulation 

Fig. 5. — The solid black line shows the projection of the 2- 
dimensional probability distribution function (the gray-scale of 
Figure [4{ projected onto the e-axis. The dashed line shows the 
e-distribution if one requires that all planets have Bond albedos 
less than 0.1; under this assumption, we see hints of a bimodal 
distribution in heat circulation efficiency. 




Albedo 

Fig. 6. — The solid black line shows the projection of the 2- 
dimensional probability distribution function (the gray-scale of Fig- 
ure |4)l projected onto the A^-axis. The cumulative distribution 
function (not shown) yields a lrr upper limit of Ag < 0.35. 

The solid line in Figure [5] shows no evidence of bi- 
modality in heat redistribution efficiency, although there 
is a wide range of behaviors. The dashed line in Figure[5] 
shows the e-distribution if one requires the albedo to be 
low, Ab < 0.1. There are then many high-recirculation 
planets, since advection is the only way to depress the 
day-side temperature in the absence of albedo. Inter- 
estingly, the dashed line does show tentative evidence of 
two separate peaks in e: if short-period giant planets 
have uniformly low albedos, then there appear to be two 
modes of heat recirculation efficiency We revisit this 
idea below. 

Figure [6] shows that planets in this sample are consis- 
tent with a low Bond albedo. Note that this constraint 
is based entirely on near- and mid-infrared observations, 
and is thus independent from the cla ims of low albedo 
based on searches for reflected light (|Rowe et al.l I2008L 
and references therein). Furthermore, this is a constraint 
on the Bond albedo, rather than the albedo in any lim- 
ited wavelength range. 

In Figure [7] we plot the dimensionless day-side effec- 
tive temperature, T^/Tq, against the maximum expected 
day-side temperature, T £= q. The cyan asterisks in Fig- 
ure [7] show the four hot Jupiters without temperature 
inversion s, while most of the remaining planets have in- 
versions (|Knutson et al.ll2010D . The presence or absence 
of an inversion does not appear to affect the efficiency of 
day-night heat recirculation. 

Planets should lie below the solid red line in Figure [3 
which denotes T £=0 = (2/3) 1 / 4 T . Of the 24 planets in 
our sample, only one (Gl 436b) has a day-side effective 
temperature significantly above the T e= q limifl This 
planet is by far the coolest in our sample, it is on an ec- 
centric orbit, and observations i ndicate that it may hav e 
a non-equilibrium atmosphere (St evenson et all 120101 ). 
There is no reason, on the other hand, that planets 
shouldn't lie below the red dotted line in Figure [7J 
all it would take is non-zero Bond albedo. That said, 
only 3 of the 24 planets we consider are in this region, 

6 This is driven by the abnormally high 3.6 micron brightness 
temperature; including the 4.5 micron eclipse upper limit does not 
significantly change our estimate of this planet's effective temper- 
ature. 
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Fig. 7. — The dimensionless day-side effective temperature, 
Td/To, plotted against the maximum expected day-side temper- 
ature, T e= o. The red lines correspond to three fiducial limits of 
recirculation, assuming Ag = 0: no recirculation (solid), uniform 
day- hemisphere (dashed), and uniform planet (dotted). The gray 
points indicate the default values (using only observations with 
A > 0.8 micron) for the four planets whose optical eclipse depths 
may be probing thermal emission rather than just reflected light 
(from left to right: TrES-2b, CoRoT-2b, CoRoT-lb, HAT-P-7b). 
For these planets we have here elected to include optical measure- 
ments in our estimate of the day-side bolometric flux and effective 
temperature, shown in black. The cyan asterisks denote those hot 
Ju piters known not to have a stratospheric inversion according 
to IIKnutson et al J 120101) . They are, from left to right: TrES-lb, 
HD 189733b, TrES-3b, WASP-4b. The two red x's denote the ec- 
centric planets in our sample, which are also the two worst outliers. 

with the greatest outlier being HD 80606b, a planet on 
an extremely eccentric orbit with superior conjunction 
nearly coinciding with periastron. As such, it is likely 
that much of the energy absorbed by the planet at that 
point in its orbit performs mechanical wo rk (speeding up 
wind s, puffing up the planet, etc. See also lCowan fc Agol 
2010) rather than merely warming the gas. Gl 436b and 
HD 80606b are denoted by red x's in Figure [7J 

The gray points in Figure [7] indicate the default val- 
ues (using only observations with A > 0.8 micron) for the 
four planets whose optical eclipse depths may be probing 
thermal emission rather than just reflected light (from 
left to right: TrES-2b, CoRoT-2b, CoRoT-lb, HAT- 
P-7b). For these planets we have here elected to use 
all available flux ratios (including optical observations 
potentially contaminated by reflected light) to estimate 
the day-side bolometric flux and effective temperature, 
shown as black points in Figure [7] 

If one takes these day-side effective temperature es- 
timates at face value, it appears that the planets with 
T £= q < 2400 K exhibit a wide- variety of redistribution 
efficiencies and/or Bond albedos, but are consistent with 
Ab = 0. It is worth noting that many of the best char- 
acterized planets in this region have T^/Tq ss 0.75, and 
this accounts for the sharp peak in the dotted line of Fig- 
ure O at e — 0.75. The hottest 6 planets, on the other 
hand, have uniformly high T^/Tq, indicating that they 
have both low Bond albedo and low redistribution effi- 
ciency. These planets must not have th e high-altitude, 
reflect ive silicate clouds hypothesized in ISudarskv et al.l 
(2000). But this conclusion is dependent on how one 
interprets the Kepler observations of HAT-P-7b: if the 
large optical flux ratio is due to reflected light, then this 
planet is cooler than we think, and even the hottest tran- 



siting planets exhibit a variety of behaviors. 

5. SUMMARY & CONCLUSIONS 

We have described how to estimate a planet's incident 
power budget (To), where the uncertainties are driven by 
the uncertainties in the host star's effective temperature 
and size, as well as the planet's orbit. We then described 
a model-independent technique to estimate the effective 
temperature of a planet based on planet/star flux ra- 
tios obtained at various wavelengths. When the observed 
day-side and night-side effective temperatures are com- 
pared, one can constrain a combination of the planet's 
Bond albedo, Ab, and its recirculation efficiency, e. We 
applied this analysis on 24 known transiting planets with 
measured infrared eclipse depths. 

Our principal results are: 

1. Essentially all of the planets are consistent with low 
Bond albedo. 

2. We firmly rule out the "null hypothesis", whereby all 
transiting planets can be fit by a single Ab and e. It 
is not immediately clear whether this stems from differ- 
ences in Bond albedo, recirculation efficiency, or both. 

3. In the few cases where it is possible to unambiguously 
infer an albedo based on optical eclipse depths, they are 
extremely low, implying correspondingly low Bond albe- 
dos (< 10%). If one adopts such low albedos for all 
the planets in our sample, the discrepancies in day-side 
effective temperature must be due to differences in recir- 
culation efficiency. 

4. These differences in recirculation efficiency do not 
appear to be correlated with the presence or absence of 
a stratospheric inversion. 

5. Planets cooler than T e=u = 2400 K exhibit a wide va- 
riety of circulation efficiencies that do not appear to be 
correlated with equilibrium temperature. Alternatively, 
these planets may have different (but generally low) albe- 
dos. Planets hotter than T e= o = 2400 K have uniformly 
low redistribution efficiencies and albedos. 

The apparent decrease in advective efficiency with 
increasing planetary temperature remains unexplained. 
One hypothesis, mentioned earlier, is that TiO and VO 
would provide additional optical opacity in atmospheres 
hotter than T ~ 1700 K, leading to temperature in- 
vers ions and reduced he at recirculation on these plan- 
ets (jFortnev et al.ll2~008h . But if our sample shows any 
sharp change it behavior it occurs near 2400 K, rather 
than 1700 K. One could invoke another optical absorber, 
but in any case the lack of correlation - pointed out in 
this work and elsewhere — between the presence of a tem- 
perature inversion and the efficiency of heat recirculation 
makes this explanation suspect. Another possible expla- 
nation for the observed trend is that the hottest planets 
have the most ionized atm ospheres and may suffer the 
most severe magnetic drag (|Perna et al.ll2010l ). 

The simplest explanation for this trend is simply that 
the radiative time is a steeper function of temperature 
than the advective time: advective efficiency is given 
rough ly by the ratio of the radiative and advective times 
(eg: iCowan &; Ago] I2010T) . In the limit of Newtonian 
cooling, the radiative time scales as r rac i a T" 3 . If one 
assumes the wind speed to be of order the local sound 
speed, then the advective time scales as r a d v oc T~°- a . 
One might therefore naively expect the advective effi- 
ciency to scale as T . Such an explanation would not 
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explain the apparent sharp transition seen at 2400 K, 
however. 

The combination of optical observations of secondary 
eclipses and thermal observations of phase variations is 
the best way to constrain planetary albedo and circu- 
lation. The optical observations should be taken near 
the star's blackbody peak, both to maximize signal-to- 
noise, and to avoid contamination from the planet's ther- 
mal emission, but this separation may not be possible for 
the hottest transiting planets. The thermal observations, 
likewise, should be near the planet's blackbody peak to 
better constrain its bolometric flux. Note that this wave- 
length is shortward of the ideal contrast ratio, which typ- 
ically falls on the planet's Rayleigh- Jeans tail. Further- 
more, the thermal phase observations should span a full 
planetary orbit: the light curve minimum is the most sen- 
sitive measure of s, and should occur nearly half an orbit 
apart from the light cu rve maximum, despite skew ed di- 
urnal heating patterns ()Cowan AAgpl 2008, 2010). This 
means that observing campaigns that only cover a little 
more than half an orbit (transit — » eclipse) are probably 
underestimating the real peak-trough phase amplitude. 



A possible improvement to this study would be to per- 
form a uniform data reduction for all the Spitzer exo- 
planet observations of hot Jupiters. These data make up 
the majority of the constraints presented in our study 
and most are publicly available. And while the pub- 
lished observations were analyzed in disparate ways, a 
consensus approach to correcting detector systematics is 
beginning to emerge. 
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TABLE 1 

Secondary Eclipses & Phase Variations of Exoplanets 



Planet 


T £=0 [K] a 


A [fim] b 


Eclipse Depth c 


^bright [K] 


Phase Amplitude 2 


Derived Quantities' 1 


CoRoT-lb 1 


2424(84) 


0.60(0.42) 


1.6(6) x 10~ 4 


2726(141) 






T d =2674(144) K 






0.71(0.25) 


1.26(33) x 10~ 4 


2409(75) 


1.0(3) X 10~ 4 




A\ <0.1 






2.10(0.02) 


2.8(5) x 10~ 3 


2741(125) 




T d (A 


= 0)=2515(84) K 






2.15(0.32) 


3.36(42) x 10~ 3 


2490(157) 












3.6(0.75) 


4.15(42) x 10~ 3 


2098(116) 












4.5(1.0) 


4.82(42) x 10~ 3 


2084(106) 








CoRoT-2b^ 


1964(42) 


0.60(0.42) 


6(2) x 10~ b 


2315(85) 






T d =1864(233) K 






71 (0 25^1 


1.02(20) x 10~ 4 


2215f49 = ) 






A\ = 0.16(7) 






1.65(0.25) 


< 1.7 x 10~ 3 (3<r) 






T d (A = 


= 0)=2010(144) K 






2.15(0.32) 


1.6(9) x 10~ 3 


1914(292) 












3.6(0.75) 


3.55(20) x 10~ 3 


1798(40) 












4.5(1.0) e 


4.75(19) x 10~ 3 


1791(33) 












4.5(1.0) 


5.10(42) x 10~ 3 














8.0(2.9) 


4.1(1.1) x 10~ 3 














8.0(2.9) e 


4.09(80) x 10~ 3 


1318(143) 








Gl 436b' 5 


934(41) 


3.6(0.75) 


4.1(3) x 10~ 4 


1145(23) 






T d =1082(38) K 






4.5(1.0) 


< 1.0 X 10~ 4 (3a) 














5.8(1.4) 


3.3(1.4) X 10~ 4 


797(106) 












8.0(2.9) e 


4.52(27) x 10~ 4 


737(17) 












8.0(2.9) 


5.7(8) x 10~ 4 














8.0(2.9) 


5.4(7) x 10~ 4 














16(5) 


1.40(27) x 10~ 3 


963(126) 












24(9) 


1.75(41) x 10~ 3 


1016(182) 








HAT-P-lb 4 


1666(38) 


3.6(0.75) 


8.0(8) x 10~ 4 


1420(47) 






T d =1439(59) K 






4.5(1.0) 


1.35(22) x 10~ 3 


1507(100) 












5.8(1.4) 


2.03(31) x 10~ 3 


1626(128) 












8.0(2.9) 


2.38(40) x 10" 3 


1564(151) 








HAT-P-7b b 


2943(95) 


0.65(0.4) 


1.30(11) x 10 -4 


3037(35) 


1.22(16) x 10 -4 




T d =2086(156) K 






3.6(0.75) 


9.8(1.7) x 10~ 4 


2063(152) 






A\ = 0.58(5) 






4.5(1.0) 


1.59(22) x 10~ 3 


2378(179) 




T d (A 


= 0)=2830(86) K 






5.8(1.4) 


2.45(31) x 10" 3 


2851(235) 












8.0(2.9) 


2.25(52) x 10~ 3 


2512(403) 








HD 80606b° 


1799(50) 


8.0(2.9) 


1.36(18) x 10~ 3 


1137(73) 




T d =1137(113) K 


HD 149026b' 


1871(17) 


8.0(2.9) e 


3.7(0.8) x 10~ 4 


976(276) 


2.3(7) x 10~ 4 




T d =1571(231) K 






8.0(2.9) 


8.4(1.1) x 10~ 4 








T n =976(286) K 


HD 189733b B 


1537(16) 


2.15(32) 


< 4.0 x 10~ 4 (1<t) 








T d =1605(52) K 






3.6(0.75) 


2.56(14) x 10~ 3 


1639(34) 






T n =1107(132) K 






4.5(1.0) 


2.14(20) x 10~ 3 


1318(45) 












5.8(1.4) 


3.10(34) x 10~ 3 


1368(69) 












8.0(2.9) 


3.381(55) x 10~ 3 














8.0(2.9) 


3.91(22) x 10~ 3 




1.2(2) x 10~ 3 










8.0(2.9) e 


3.440(36) x 10~ 3 


1259(7) 


1.2(4) x 10~ 3 












5.51(30) x 10~ 3 


1338(52) 












24(9) 


5.98(38) x 10~ 3 














24(9) e 


5.36(27) x 10~ 3 


1202(46) 


1.3(3) x 10" 3 






HD 209458b' J 


1754(15) 


0.5(0.3) 


7(9) x 10 _t> 


2368(156) 






T d =1486(53) K 






2.15(0.32) 


< 3 x 10- 4 (1(t) 








A\ = 0.09(7) 






3.6(0.75) 


9.4(9) x 10~ 4 


1446(45) 




T d (A = 


= 0)=2031(128) K 






4.5(1.0) 


2.13(15) x 10~ 3 


1757(57) 






T n =1476(304) K 






5.8(1.4) 


3.01(43) x lO" 3 


1890(149) 












8.0(2.9) 


2.40(26) x 10~ 3 


1480(94) 


< 1.5 x 10~ 3 (2a) 










24(9) 


2.60(44) x 10~ 3 


1131(143) 








OGLE-TR-56b iu 


2874(84) 


0.90(0.15) 


3.63(91) x 10~ 4 


2696(116) 




T d =2696(236) K 


OGLE-TR-llSb 11 


1716(33) 


2.15(0.32) 


1.7(5) x 10~ 3 


1918(164) 




T d =1918(219) K 


TrES-lb^ 


1464(16) 


3.6(0.75) 


< 1.5 x 10~ 3 (1<t) 








T d =998(67) K 






4.5(1.0) 


6.6(1.3) x 10~ 4 


972(56) 












8.0(2.9) 


2.25(36) x 10~ 3 


1152(94) 








TrES-2b iJ 


1917(21) 


0.65(0.4) 


1.14(78) x 10 _0 


2020(132) 






T d =1623(76) K 






2.15(0.32) 


6.2(1.2) x 10~ 4 


1655(80) 






A\ = 0.06(3) 






3.6(0.75) 


1.27(21) x 10~ 3 


1490(84) 




T d (A = 


= 0) = 1751(80) K 






4.5(1.0) 


2.30(24) x 10~ 3 


1652(74) 












5.8(1.4) 


1.99(54) x 10" 3 


1373(177) 












8.0(2.9) 


3.59(60) x 10~ 3 


1659(163) 








TrES-3b 14 


2093(32) 


0.7(0.3) 


< 6.2 x 10~ 4 (la) 








T d =1761(66) K 






1.25(0.16) 


< 5.1 X 10~ 4 (3<rj 














2.15(0.32) 


2.41(43) X 10~ 3 














2.15(0.32) e 


1.33(17) x lO" 3 


1770(58) 












3.6(0.75) 


3.46(35) x 10" 3 


1818(73) 
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4.5(1.0) 


3.72(54) 


x 10~ 3 


1 649 ("1 071 








5.8(1.4) 


4.49(97) 


x 10~ 3 


1621(173) 








8.0(2.9) 


4.75(46) 


x 10~ 3 


1480(82) 




TrES-4b lb 


2250f37l 




1.37(11) 


x 10~ 3 


1 889("fi3) 

iou3i yj%j j 


Tj — 18Q1 (81 1 K 

J- lUC/ll Ol 1 IV 






4.5(1.0) 


1.48(16) 


x 10~ 3 


1 727(831 








5.8(1.4) 


2.61(59) 


x 10~ 3 


2112(283) 








8.0(2.9) 


3.18(44) 


x 10~ 3 


2168(197) 




WASP-lb ib 


2347fS5l 


3 6(0 75"! 


1.17(16) 


x 10~ 3 


1678(87) 


Ti— 1719(891 K 

J. ^ — i j ici i iv 






4.5(1.0) 


2.12(21) 


x 10~ 3 


1923(91) 








5.8(1.4) 


2.82(60) 


x 10~ 3 


2042(253) 








8.0(2.9) 


4.70(46) 


x 10~ 3 


2587(176) 






1 fifil f69l 


3 fiCO 75") 


8.3(3.5) 


x 10~ 4 


1 264(1 64"! 


Tj— 1280(1211 K 

_1 {-] lijUUl 1^.1 1 IV 






4.5(1.0) 


1.69(17) 


x 10~ 3 


1 380(531 








5.8(1.4) 


1.92(77) 


x 10~ 3 


1 299(2321 








8.0(2.9) 


2 85(591 


x 10~ 3 


1 372(1 541 




WASP 4h itt 


91 fi'3ffifj^ 


^ fii'n 7^ ) 

O.U y\J . 1 O ) 


3.19(31) 


x 10 — 3 


91 e: if\(Q7\ 

£ -LvJU \ J t J 


± J — £ ±<-k\jy ±^±U ^ IV 






4.5(1.0) 


3.43(27) 


x 10~ 3 


1 Q71 ^7^ 




WA SP 1 9K iy 

V V Jr\ Or - 1ZL) 


ozio^iiy ) 




o. z ^±.U^ 


x 10 -4 


Qnfl9f 1 f\A\ 


T 1 ! — 9Q'3Q/'QR,' K 
id — zyoy^yo j rv 








l .Ol J 


X 10 — 3 










i fs^?n 9^ , 


1.76(18) 


x 10~ 3 


9«9Q('e«^ 

iiOinJ luO / 








2.15(0.32) 


3.09(13) 


x 10~ 3 


3018(51) 








3.6(0.75) 


3.79(13) 


x 10~ 3 


2704(49) 








4.5(1.0) 


3.82(19) 


x 10~ 3 


2486(68) 








5.8(1.4) 


6.29(52) 


x 10~ 3 


3167(179) 








8.0(2.9) 


6.36(67) 


x 10~ 3 


2996(229) 




WASP-18b^ u 


3070(50) 


3.6(0.75) 


3.1(2) 


x 10~ 3 


3000(107) 


T d =2998(138) K 






4.5(1.0) 


3.8(3) 


x 10~ 3 


31280501 








5.8(1.4) 


4.1(2) 


x 10~ 3 


3095(103) 








8.0(2.9) 


4.3(3) 


x 10~ 3 


2991(153) 




WASP-lQb^ 1 


2581(49) 


1.65(0.25) 


2.59(45) 


x 10 -3 


2677(135) 


T d =2677(244) K 


XO-lb-" 


1 t S2fi('24'l 

J. tJZj VJ 1 ^'i 1 


3 6("0 75") 


8.6(7) 


x 10 -4 


1 300(321 


Tj-n0flf47, K 

J. d -LOVJVJtril J IV 






4.5(1.0) 


1.22(9) 


x 10~ 3 


1 265 ( 341 








5.8(1.4) 


2.61(31) 


x 10~ 3 


1546(89) 








8.0(2.9) 


2.10(29) 


x 10" 3 


1211(87) 






1685(33) 


3.6(0.75) 


8.1(1.7) 


x 10- 4 


1447(102) 


T d =1431(98) K 






4.5(1.0) 


9.8(2.0) 


x 10~ 4 


1341(105) 








5.8(1.4) 


1.67(36) 


x 10~ 3 


1497(155) 








8.0(2.9) 


1.33(49) 


x 10" 3 


1179(219) 




XO-3 M 


1982(82) 


3.6(0.75) 


1.01(4) 


x 10" 3 


1875(30) 


T d =1871(63) K 






4.5(1.0) 


1.43(6) 


x 10~ 3 


1965(40) 








5.8(1.4) 


1.34(49) 


x 10~ 3 


1716(330) 








8.0(2.9) 


1.50(36) 


x 10~ 3 


1625(236) 





The planet's expected day-side effective temperature in the absence of reflection or recirculation (Ab — 0, £ — 0). The la uncertainty is shown 
in parcnthcsc. 

k The bandwidth is shown in parenthese. 

C Eclipse depths and phase amplitudes arc unitlcss, since they are measured relative to stellar flux. 

^ Td and T n denote the day-side and night-side effective temperatures of the planet, as estimated from thermal secondary eclipse depths and 
thermal phase variations, respectively. The estimated la uncertainties arc shown in parentheses. The default day-side temperature is computed 
using only observations at A > 0.8 fim. Eclipse measurements at shorter wavelengths may then be used to estimate the planet's albedo at those 
wavelengths, A\. Note that this is a spherical albedo; the geometric albedo is given by A g — ^A\. If — on the other hand — Ab — is assumed, 
then all the day-side flux is thermal, regardless of waveband, yielding the second Td estimate. 

6 When multiple measurements of an eclipse depth have been published in a given waveband, we use the most recent observation. In all cases 
these observations are either explicitly agree with their older counterpart, or agree with the re-analyzed older data. 

1 ISnellen et~aT1 | |2009| '); [ATonso et all < |2009by) ; [Gillon et all 1 120091) ; [Rogers et al.| j2009D ; [Deming et all J2010f) , lAlonso et~aT| | |2009af) ; [Snellen et aT| 
l|2010D : IGillon et all d2010f ): lAlonso et~aTI d2010f) : iDeming et al.1 i2010l) . iDeming et al.l i2007D : iDemorv et all d2007D : IStevenson et all d201d 1: 
Knutson et al. in prep., iTodorov et all d2010f) . iBorucki et al.l i2009l ): IChristiansen et all i2010f) . * lLaughlin et alj i2009l) . iKnutson et alj 
d2009b!) . iDeming et all i2006f) : IKnutson et all d2007af) : I Barnes et"aTI i2007fl : ICharbonneau et al.1 J2008I 1: IKnutson et alj i2009cl) : lAgol et all 
J2010ll , c Richardso n et alj j2003f ) ; [Deming et al.1 j2005l) : [Cowari et al.l 1120070 : [Rowe et alj J2001) ; [Knutson et all i200S|) . lc ISing Lopez-Moralea 
1 120091) . 1 ^Snellen fc Covinol 120071) . ICharbonneau et alj ( 120051) : IKnutson eFal] d2007bl) . HO'Donovan et alj ( 1201011 : ICroll et all I2010al ): 
IKipping &: Bakoj j |2010bl) . 1 1Fressin et al.1 ( 120101) : ICroll et al.1 ( I2010bl) : IChristiansen et al] j |2010bl) . l 1Knutson et alj j |2009al) . 16,1 ' IWheatlev et alj 
(I2010T) . ^Beerer et al.1 11201011 . ls IL6pez-Morales et aD ll2010l ): ICampo et all d2010D : rCroll et al.l<2010cl) . 2 HNvmever et al.1 d2010l) . 21 IAnderson et al.1 
(I2010T) . 2i IMachalek et al.1 <200g|) . 2 1Machalek et al.1 (120091) . 2 IMachalek et al.1 (120101 ) 



